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ABSTRACT 

The first and third year data releases from the Wilkinson Microwave Anisotropy Probe (WMAP) 
provide evidence of an anomalous Cold Spot (CS) at galactic latitude b = —57° and longitude / = 209°. 
We have examined the properties of the CS in some detail in order to assess its cosmological significance. 
We have performed a cluster analysis of the local extrema in the CMB signal to show that the CS is 
actually associated with a large group of extrema rather than just one. In the light of this we have 
re-examined the properties of the WMAP Internal Linear Combination (ILC) and co-added "cleaned" 
WCM maps, which have previously been used for the analysis of the properties of the signal in the vicinity 
of the CS. These two maps have remarkably similar properties on equal latitude rings for \b\ > 30°, as 
well as in the vicinity of the CS. We have also checked the idea that the CMB signal has a non- Gaussian 
tail, localized in the low multipole components of the signal by representing the CMB signal S outside 
the Galactic mask as a collection of signals for each equal latitude ring S = U '^(^)- ^^^r each ring we 
apply a linear filter with characteristic scale dividing the CMB signal in two parts: the filtered part, 
with characteristic scale above that of the filter and the difference between the initial and filtered 
signal. Using the filter scale as a variable, we can maximize the skewness and kurtosis of the smoothed 
signal and minimize these statistics for the difference between initial and filtered signal. We find that, 
unlike its Northern counterpart, the Southern Galactic hemisphere of the CMB map is characterized by 
significant departure from Gaussianity of which the CS is not the only manifestation: we have located a 
ring, on which there are "cold" as "hot" spots with almost the same properties as the CS. Exploiting the 
similarity of the WCM and the ILC maps, and using the latter as a guide map, we have discovered that 
the shape of the CS is formed primarily by the components of the CMB signal represented by multipoles 
between 10 < ^ < 20, with a corresponding angular scale about 5 — 10°. This signal leads to modulation 
of the whole CMB sky, clearly seen at \b\ > 30° in both the ILC and WCM maps, rather than a single 
localized feature. After subtraction of this modulation, the remaining part of the CMB signal appears 
to be consistent with statistical homogeneity and Gaussianity. We therefore infer that the mystery of 
the WMAP CS refiects directly the peculiarities of the low-multipole tail of the CMB signal, rather than 
a single local (isolated) defect or manifestation of a globally anisotropic cosmology. 



1. INTRODUCTION 

An extremely Cold Spot (CS), apparently inconsis- 
tent with the assumption of statistically homogeneous 
Gaussian fiuctuations, was detected in a wavelet analy- 
sis (Vielva et al. 2004; Cruz et al. 2005; Cayon, Jin 
and Treaster 2005; Cruz et al. 2006) of the first-year 
data release from the Wilkinson Microwave Anisotropy 
Probe (WMAP). More recently, the existence of this spot 
has been confirmed by Cruz et al. (2007a) the WMAP 
third year data release (Hinshaw et al. 2007; Spergel et 
al. 2007). The WMAP CS is centered at the position 
b = —57°, / = 209° in Galactic Coordinates and has a char- 
acteristic scale about 10°. As was pointed out by Cruz et 
al. (2006), the frequency dependence of the signal in the 
spot area is extremely flat. This fact has been used by 
the authors mentioned above to argue that the WMAP 



CS belongs to the CMB signal, rather than any form of 
foreground emission. Cruz et al. (2006) pointed out that 
the reason the CS was not been detected in real space be- 
fore the wavelet analysis was that it was hidden amongst 
structures at different scales. 

As an apparent example of non- Gaussian behaviour in 
the WMAP CMB signal, the CS has attracted very seri- 
ous attention from the theoretical point of view. Tomita 
(2005) suggested that the CS can be related to second- 
order gravitational effects. Inoue and Silk (2007) proposed 
a model involving local compensated voids. The origin 
of the CS in connection to the brightness and number of 
counts of the NVSS sources (smoothed on the scale of a 
few degrees) was recently discussed by Rudnick, Brown 
and Wilhams (2007). They have detected a 20-45% dip 
in the smoothed NVSS source counts which can be inter- 
preted, they argue, as a manifestation of the integrated 
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Sachs- Wolfe effect, seen for a single region of the CMB 
sky. Jaffe et al. (2005, 2006), Cayon et al. (2006), and 
McEwen et al. (2005, 2006) have investigated the Bianchi 
VII/i anisotropic cosmological model as a possible expla- 
nation of the CS and other features of the WMAP low 
multipoles. Recently, Cruz et al. (2007b) have pointed 
out that the CS could be produced by a cosmic texture, 
assuming that the CMB signal is a combination of the 
Gaussian and non-Gaussian parts. The present status of 
the problem of existence of the CS remains uncertain, de- 
spite the presence of a vast collection of theoretical sug- 
gestions. Nevertheless, if we believe that one particular 
part of the WMAP CMB signal contains non-Gaussian fea- 
tures, it would be necessary to seek corroborating evidence 
of non-Gaussianity elsewhere in order to understand their 
properties more fully. In this Paper we therefore present a 
detailed investigation of the properties of the CS, focusing 
attention on the following topics. 

First, in Section 2, we show how the CS can be eas- 
ily detected in the pixels domain not only in the derived 
CMB signal, but even in the WMAP maps for K-W bands 
before separation of the signal into CMB and foreground 
components. 

Second, we will demonstrate that the CS belongs to a 
cluster of local minima, the spatial distribution of which 
is modulated by the large-angle modes of the CMB signal 
outside the Galactic plane. For that we use the Inter- 
nal Linear Combination (ILC) III map and the co-added 
WCM map with Nside = 512 in the HEALPix format, 
converted to GLESP format (Doroshkevich et al. 2003), 
where each iso-latitude ring has the same number of pix- 
els N(j) = 2048 in azimuthal direction (p in polar coordi- 
nates. After that we perform a cluster analysis (Novikov 
and Jorgensen 1996) of the positive and negative peaks for 
selected rings in the area outside the KpO mask, mostly 
concentrating our attention on the ring crossing the CS at 
its extremum b = -57° and -180° < (j) < 180°. Taking 
into consideration the signal for each ring with the latitude 
6, we can investigate the morphology of the CMB signal 
at each latitude for the whole range of (j). This approach 
allows to connect the morphology of the CS to the signal 
outside the CS for the same latitude b = —57°. We will 
show that the cluster contained the CS, is not a unique fea- 
ture of the b = —57° iso-latitude ring. For example, close 
to the CS there are two significant clusters of maxima, but 
these peaks have lower amplitude that the CS. 

Next, since the origin of large clusters of extrema is 
related to the angular modulation of the signal on large 
scales (see Novikov and Jorgensen 1996), we split the CMB 
signal into two parts. To do that we use the skewness 
and the kurtosis of the signal for selected rings, includ- 
ing the b = —57° ring. Then by using a simple lin- 
ear smoothing filter with characteristic scale R we sep- 
arate the signal into a smoothed component and to a 
difference between initial signal and the smoothed com- 
ponent. For the smoothed signal we define the skewness 
S{R) and the kurtosis K{R) as a functions of R and find 
the scale of filtering which maximize both these charac- 
teristics S (Ropt) ^ K (Ropt) max. By using this scale 
Rmax we separate the initial CMB signal in two parts, 
one of them (the smoothed one) contain a maximally non- 
Gaussian signal, and the another one (initial signal minus 



the smoothed one) the maximally Gaussian signal. The 
non- Gaussian part is mainly formed by the signal local- 
ized at the range of multipoles 2 < ^ < 20 and the other 
one belongs to the ^ > 20 multipoles. Our analysis clearly 
demonstrates that the pronounced non-Gaussianity of the 
CS refiects directly the existence of a large-scale angular 
modulation of the CMB signal with 10 < ^ < 20. 

Finally, using cluster analysis in combination with skew- 
ness and kurtosis statistics, we are able to detect a few 
additional cold and hot spots on the same b = —57° ring 
as the famous one. To show that the effect on cluster- 
ing of the peak by low multipoles of the CMB is very 
common, we took into consideration the north Galactic 
hemisphere with deficit of the power and have found a few 
cold and hot spots. The idea of implementing of cluster 
analysis in combination with skewness and kurtosis statis- 
tics was stimulated by Vielva et al. (2004), Cruz et al. 
(2005, 2006, 2007a), Cayon, Jin and Treaster (2005), and 
especially Eriksen et al. (2004) and Larson and Wandelt 
(2005). 

In our analysis we use both the WCM map and the 
ILC third year map, which are very similar outside the 
b = ±25° cut of the Galactic plane. For our analysis we 
use the high resolution ILC III map as a guide map to 
mark possible zones of the CMB sky in which the enhanced 
clustering of the peaks is expected to be considerable. 



2. "naive" detection of the WMAP COLD SPOT 




Fig. 1. — The map for negative (top) and positive thresholds 
-0.2 < T < -0.05 and 0.05 < T < 0.2 of the ILC III map. For 
the top map the black circle marks the location of the CS in the 
Galactic coordinates. For the bottom map the black circles marks 
high amplitude positive peaks. 

As was pointed out in the Introduction, historically the 
CS was detected in the WMAP data as one of the deep- 
est minima of the CMB signal, using method based on 
wavelets. Our first aim is to show that this CS can actually 
be detected quite straightforwardly in the pixel domain us- 
ing simple thresholding techniques. In Fig.l we take two 
thresholds of the ILC III signal at the range of temper- 
atures -0.2 < AT < -0.05 and 0.05 < AT < O.llmK 
and map them with a color scale —0.2,0 and 0,0.2, re- 
spectively. 
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Fig. 2. The SKY Viewer maps for the ILC III (top), the K and 
Ka bands of the WMAP (second and third from the top) map with 
cutoff of intensity down to the minimal values. For all the maps the 
zone of the CS has been placed at the center of the black circle. 

From this figure one can immediately see the difference 
in statistical properties of the positive and negative thresh- 
olds of the map. First of all, it is important to note that 
a majority of the positive peaks are associated with the 
area around \b\ < 25 — 30° (see the bottom plot), while 
the negative peaks are present well outside this region of 
the map. Secondly, the amplitude of the highest positive 
peaks is limited by O.llmK, while for the negative peaks 
it is about —0.019mK . For the CS the temperature of the 
negative peak is in agreement with the estimate of Cruz et 
al. (2007a). In Fig. 2, by using the SKY Viewer, we map 
the ILC III K and Ka maps, in which the CS is located 
in the center of the area marked by the black circle (the 
blue cluster of peaks inside the black circle). From the K 
and Ka maps one can clearly see that the CS is clearly 
observed even without subtraction of the CMB signal. In 
agreement with Rudnick, Brown and Williams (2007) we 



note that there is a cluster of negative peaks, rather than 
one single peak. 

To show the local structure of the zone containing the 
CS, in Fig.3 we plot high resolution {i-max — 100) images 
of the inner structure of the CS zone, including the Haslam 
et al. (1982) map. 





Fig. 3. — The CS zone image projection (in the pixel domain) of 
the ILC III, K, Ka and Haslam et al. (1982) maps. The size of all 
the maps is 25 x 25°. 

From Fig. 2 and Fig.3 one can draw the important con- 
clusion that the zone of the CS is surrounded by zones of 
hot spots, clearly seen in Fig. 2 just on the right and left 
hand side of the CS. 

If the origin of the CS is related to large-angular mod- 
ulation of the CMB map and possible anisotropy of the 
power distribution across the sky, it would be naturally 
to expect that additionally to the WMAP CS detected by 
Vielva et al. (2004) and others, we could find other cold 
and hot spots similar or even equal morphological proper- 
ties. To show that this is the case for WCM and the ILC 
III maps, in the next section we will use a cluster analysis 
of these maps in order to mark possible zones of peculiar 
distribution (clusters) in the signal. 

3. CLUSTER ANALYSIS OF THE CMB MAPS IN THE 
VICINITY OF THE CS 

The main idea behind our implementation of cluster 
analysis is to divide the CMB sky above and below the 
KpO mask into iso-latitude rings and then to analyze the 
properties of the signal for each ring separately. First, let 
us briefiy describe the statistical properties of the signal 
T{6 = Oc,(p) for fixed latitude b. Let us take under con- 
sideration the distribution of peaks above and below some 
given threshold Ut = AT/ctq, where <Jo is the square root 
of the variance ctq for each ring. For a one-dimensional 
cross-section of the CMB map we introduce a definition of 
a cluster of maxima considering, for example, two points, 
Xi and X2 {xi < X2) from the ring. If for all the points 
in the interval xi < x < X2 we have AT{x) > z^tcro, we 
call these points xi^X2 connected to each other with re- 
spect to the threshold Uf A collection of maxima of AT 
located at the points {x^}, k = 1^2..K we will call a clus- 
ter of length D if all the points {xk} are connected to each 
other when the threshold Ut is applied ^. The length of the 
cluster D = \xi — X2|/27r is an analogue of the definition 
of the two-dimensional area (Cruz et al. 2007a), but in 
a one-dimensional case. For a random Gaussian field the 



The same definition can be introduced for cluster of minima, when instead of maxima we will use the minima of T. 
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statistic of length of the clusters (one-dimensional) and 
the area (two-dimensional) are similarly sensitive to the 
spectral parameters of the random field, as described by 
by Bardeen et al. (1986). However, the diff'erence between 
cluster analysis and area statistics becomes obvious in ap- 
plication to non- Gaussian random fields. Let us assume 
for a moment that some zone of the signal is characterized 
by the area statistic, and that this has a very non- Gaussian 
value. To characterize the properties of the non- Gaussian 
field it seems to be very important to know, if this area 
is related to a uniform structure with one single peak, or 
formed from a cluster of peaks with the same area. This is 
why for our analysis of the WCM and the ILC III statis- 
tical properties we prefer to use cluster analysis (CA). In 
our analysis, in addidition to the length of cluster D, we 
will use its size S which we define as the number of max- 
ima (or minima) above (or below) the threshold ut within 
the interval length D = \xi — X2\/27r. 

3.1. Statistical properties of the signals for equal latitude 

ring 

To describe the statistical properties of equal latitude 
rings we will use the approach proposed by Chiang and 
Naselsky (2007). 

The standard treatment for a full-sky CMB signal 
T{0,(p) is via spherical harmonic decomposition: 

T{0,^) = J2 E airuYemiO,^), (1) 

i=0 m=-i 

where and cp are the polar and azimuthal angle, respec- 
tively, and aim are the spherical harmonic coefficients. 
Here are the spherical harmonics, defined in terms 
of Legendre polynomials and plane waves: 

Yim{0, ^) = NimPn^osO) exp(im99), (2) 

where 



= {-ly 



l{2£^l){£-m)\ 



(3) 



For a random Gaussian CMB sky the properties of the 
statistical ensemble of realizations are completely specify 
by the power spectrum 



c{e) = 



2£-hl ^ 



ml )^ 



(4) 



while for each single realization we expect to find some 
deviation from C{£) due to the "cosmic variance" effect. 
Below we will use a polar coordinate system in which the 
Galactic plane (b = 0) is associated with = 7r/2. 

Let us analyse the signal TiOc^if) from the equal- 
latitude ring at = 0c, where \0c\ > |6>mask|, where |6>mask| 
is the maximum latitude of any foreground masks. This 
ring T(^c, ^) = ^c(^) is a one-dimensional signal, for 
which we can construct a Fourier transform with coefli- 
cients gm- 



gm exp(im(/?). 



where 



gm= dipTc{(p) exp{-imip). 
Jo 



(5) 



(6) 



We can then relate the ring to the full-sky signal via Eq.(2) 
and (5) and get 

-''max 

9m= (^^^ • ( 

i>\m\ 

That is to say that the Fourier coefficients gm of the ring 
can be expressed as a combination of the full-sky aim- 
Defining the variance of the signal for equal latitude ring 
as 

(8) 

and after substitution of Eq.(5) into Eq.(8) we have 

{T) = Ni^m=oPi{cos O)ai,m=o, 
i 

xPi^^m{cOsO)ai^mCi}r,m' (9) 

For a random Gaussian field (GRF), after average over 
realizations, the combinations of the aim coefficients in 
Eq.(9) satisfy the following conditions: 

{(^£,mCit,m') = C{£)6£/^6m,m' (10) 

and then 

VarT = J2Y1 Nl^Plm{cos0)C{£) = ^W{£,cose)C{l) 



i m/0 



(11) 



where W{£,cosO) = Y^m^o ^ImPlmi^^^^) ^^le window 
function of the ring. 

Our approach here is a special case for a well known 
theoretical prediction: any n dimensional cross sections of 
N dimensional Gaussian random signal produce a Gaus- 
sian signal as well. However, this general theory tells us 
about the GRF, represented as a sets of realizations. For 
one single realization of the CMB sky, as the WMAP sig- 
nal does, the Eq.(ll) is no longer available. More general, 
instead of Eq.(ll) we have 

«^,m^i|/,m' = ^{^)^T,Jt^ ' i^T,Jt^ ) ^ ^i,i'^m,m' (12) 

Thus, 

,m'Pe,m{cOsO) X 

Pe,^m'{cose)C{l)Glf , 

,m' X 

i' m/0 

P,,^(cos^)P,,,^,(cos^)G^^r' (13) 

The matrix G^'^ describes the coupling between differ- 
ent modes £,m and £' ^m' for the random process T((/?), 
which leads to variations of the variance VarT{0) for 
different rings {0 = const). One can see that integra- 
tion of Eq.(13) over preserves spontaneous correlations 

(G^^J^ 7^ S£/'Sm,m') even for the whole sky. However, for 
a particular ring we will have an additional modulation of 
these correlations depending on the through the window 
function l^(^,cos^). In Fig. 4 we show the dependence of 
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VarT on the galactic latitude for the WCM and the ILC 
III maps. 




-50 50 
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Fig. 4. — The variance for the WCM map (the black line) and 
the ILC III map (the red line) versus the galactic latitude. 

First of all, we stress the remarkable similarity of these 
maps outside the region of the KpO mask. Moreover, the 
variance of the WCM map per each ring matches the vari- 
ance of the ILC III almost exactly for all latitudes, except 
for the \b\ < 5° zone. Secondly, note that the b = —57° 
ring lies near to the local maxima of the variance. The 
width of this zone is about ~ 10°. Thirdly, from Fig. 4 
one can clearly see the asymmetry of the variance for the 
rings located at the North and the South hemispheres. No 
special analysis is required to observe that the variance 
of the Southern hemisphere is higher than for the North. 
This result is in agreement with Eriksen et al. (2004) inves- 
tigations of the asymmetry of the power per 10° patches of 
the sky, and Larson and Wandelt et al. (2005) analysis of 
peak statistics. Furthermore, looking at Fig. 4 at 6 = —25° 
one can find the global maxima of the variance outside the 
KpO mask. This zone is clearly seen in Fig.l just below 
the galactic plane as a large cluster of minima. Could the 
origin of this cluster be the same as for the cluster around 
the CS ? To answer this question we need to look closely 
on the properties of the signal for the ring b = —57° in the 
azimuthal direction. 



3.2. The clusters in the ring at b = —57° 

In this section we draw attention to the azimuthal distri- 
bution of the signal for the iso-latitude ring with b = —57°, 
which contains the CS. In Fig.5 we plot T(6>c, (/>) for WCM 
and the ILC III maps smoothed by the angle 0c = 1°. 
This figure clearly demonstrates that there are no signif- 
icant differences of morphology of the CS in the ILC III 
and WCM maps. In the following analysis we change the 
reference system of coordinate from the Galactic one to 
one in which ^ = is associated with the Galactic center 
and then all the values of (j) are counted clockwise up to 
(j) = 360°. The GLESP pixelization allow us to fix the same 
number of pixels for each iso-latitude ring {N = 2048) and 
consequently the location of each pixel < /c < A/" is re- 
lated to the angle ^/^ as ^/^ = 27rk/N. 



0.2 
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Fig. 5. — Plot of the one-dimensional ring through the center 
of the CS (6 = —57°) in azimuthal direction. The solid black line 
corresponds to the WCM signal and the red line is for the ILC III 
map. All the maps were smoothed by the angle Oc = 1°. For 
this ring ctq = 0.0651mK. The dashed lines mark the thresholds 
Wt\ = 1,2,3. 

In Fig.5 the CS is clearly seen at /c 1130 ^ 1224 as a 
cluster of minima (CIO, n=10) of the signal in the WCM 
and ILC III maps. Moreover, there are no significant dif- 
ferences between these two signals either for the whole 
ring, or for the zone of the CS. Fig. 6 shows the distri- 
bution of the length of the clusters versus the number of 
cluster. Clusters of maxima are shown in black, minima in 
red. As one can see from Fig. 6 there are about 17 clusters 
of maxima and 17 clusters of minima. The mean length of 
the ILC and WCM clusters of maxima {D^^{n)) 0.0337, 
while for the clusters of minima we get {D^^{n)) 0.0256. 
For comparison, after taking 1000 realizations for the GRF 
with the WMAP best fit ACDM power spectrum we have 
found (I^+f(n)) 0.0303 and {D-^{n)) 0.0315. 

The cluster CIO containing the CS has length Dqs = 
0.0459 :^ 1.5(I)~f(n)) and its dimension is D~ = 2 at 

Ut = 0. At the same time Dqs/ {D^^{n)) 1.8 at Ut = 0. 
These ratios tell us that this cluster looks like an "ordi- 
nary" cluster defined by the threshold Vt = ^ without any 
significant deviation from the statistical properties of a 
random Gaussian Field. However, the extraordinary prop- 
erties of this cluster become obvious when we take into 
consideration other thresholds, —3 < Vt ^ —2. From Fig. 6 
we see that even for this range of Vt the cluster still ap- 
pears as a cluster of minima at the level v = —3, and local 
maxima at Ut = —2. Novikov and Jorgensen (1996) have 
pointed out that for the cluster of minima with d~ = 2 
the conditional probability P(/ci, /C2) to find two negative 
peaks located at the points ki and k2 in this cluster is 
less than the conditional probability to find one maxima 
in between these points with height above In out case 
this threshold is Ut ^ -2 and P(/ci, ^2) < 0.02 

However, this cluster is not a unique feature of the 
b = —57° ring. There are two clusters of maxima (see 
Fig.5), namely the cluster C9 with kmin = 825 and 
femax = 1129 with Dc9 = 0.15 = 4.40(L>+^(n)) and 
the cluster Cll with kmin = 1310, /cmax = 1634, and 
Dcii = 0.15 = 4.69(L)^:^^(n)). Here /cmin, /^^max mark the 
coordinates of the cluster at = in the pixel domain. 
These clusters seen in Fig. 6 are near to the cluster of min- 
ima CIO within which the CS lies. 
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Fig. 6. — Top. The length D(n) of the clusters i/t = versus 
the number of cluster. The black line is for clusters of maxima, the 
red line represents the minima. Bottom plot. Number of cluster 
versus its position in the pixel domain (2048 pixels correspond to 
the 360° in azimuthal direction). The blue line corresponds to a 
one-dimensional ring of the ILC III map trough the CS in units of 
5 * Tj^/dQ. The negative sign for D~ (n) marks the length of the 
clusters of minima. 

The existence of such clusters with lengths above 
4{D^^{n)) and the cluster of minima CIO with such high 
amplitude (negative) peaks seems to be a quite peculiar 
feature of the b = —57° ring. However, these peculiarities 
are not so specific. To show that "over-clustering" of the 
CMB sky is a typical feature of the signal let us examine, 
for example, the ring with b = +57°, located symmetri- 
cally to the b = —57° ring in respect to the Galactic plane. 
In Fig. 7 we show T(/c), and the distribution of clusters 
versus their length for the ring with b = 57°. This signal 
reveals a remarkable similarity of the morphology to the 
ring b = —57°. As for the ring with CS, the ring b = 57° is 
characterized by very high level of clusterisation, the exis- 
tence of the 3(7 minima, as a member of the cluster with 
D~ = 4, and the existence of the positive cluster with 
= 6. Moreover, one can see that the cluster of max- 
ima C7 at kmin = 547, k^ax = 633 has a structure similar 
to the structure of the signal in the CS, but now for the 
maxima. 

At the end of this section we would like to point out that 
the existence of clusters with the length D > 3 + 4(1)) is a 
quite rare event for the GRF. The presence of three clus- 
ters with D ~ 3{D) for single iso-latitude ring at 6 = 57°, 
and two clusters with D > 4{D) for b = -57°, as the CS 
as well, allow us to conclude that the ILC III and WCM 
maps are generally "over-clustered" . Attention on the CS 
was focused mainly because of the amplitude of the signal 
in that position. A more specific feature of the CS is not 
this, but that it is really a cluster of peaks with nearly the 
same amplitude, and there is a very large cluster near to 
it. 
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Fig. 7. — The same as for Fig. 5 and Fig. 6, but for the ring with 
6 = 57°. 



4. ''DE-CLUSTERING" OF THE ILC III AND WCM BY 
LINEAR FILTRATION 

In this Section we re-examine an idea of Cruz et al 
(2007a), that the origin of the CS can be explained if 
the WMAP CMB signal contains non-Gaussian compo- 
nents in combination with the Gaussian ones. This idea 
seems to be quite natural since the low multipole part of 
the WMAP ILC III map reveals significant peculiarities of 
the signal: alignment between quadrupole and octupole, 
coupling with the Galactic foregrounds, low power of the 
quadrupole, and so on. Note that the WMAP team per- 
formed the analysis of Gaussianity of the CMB signal by 
subtracting the ten multipole components of lowest order 
from the map and claimed that the rest of the signal is 
in agreement with Gaussian statistics. However, asymme- 
try of the power of the CMB, discovered by Eriksen et 
al. (2004), the existence of the CS and another peculiar- 
ities, mentioned in Lasenby et al. (2007) and McEwen et 
al. (2007), raise the question of whether there is some 
mark imprint non- Gaussian features of the CMB, mainly 
localised in the low multipole range of the power spectrum. 

If the WMAP CMB outside the Galactic mask |6| = 25° 
is the sum of Gaussian and non-Gaussian components, it 
seems natural to guess that these two signals would have 
different characteristic scales. The idea, which we will 
develop below, is to use a linear filter of the CMB sig- 
nal for each ring of the map with variable scale of filtra- 
tion which can divide the CMB signal in two parts: 
S{k) = S{k^R) + s(/c), where s{k) corresponds to Gaus- 
sian component and G{k) = S{k^ R) corresponds to non- 
Gaussian one. Note that Cruz et al. (2007) use the wavelet 
approach to perform this analysis. We will instead use the 



7 



running window filter defined in the pixels domain as 

R-l 

s{k,R) = {n) = -Y,T,+j-R, 

3=0 

Sk=Tk- {Tk{R)). 

(14) 

Note that the choice of the linear filter is not so important 
for the criteria of separation of the Gaussian and non- 
Gaussian tails of the signal. One can for example use the 
Gaussian filter S{k) ex Tj exp [—{k — or any 

another reasonable filter. What is important is the chose 
of the criteria of the scale R. Following Cruz et al. (2006), 
we will use the skewness and kurtosis of the signals 5/e as 
a functions of R trying to minimize the difference between 
their actual values and most probable values for the GRF. 

In Fig. 8 we show the skewness and kurtosis for the ILC 
III ring b = —57° for different scales of filtering R. 
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Fig. 8. — Top. The skewness and kurtosis of the ILC III b = -57° 
ring. The sohd black and the dash black lines corresponds to the 
skewness and kurtosis of the smoothed part of the ILC III signal 
S(k,R). The red lines are for the signal Sk- Bottom. The same as 
top, but for the zone of the ring kmin = 825, kmax = 1634. 

The skewness and kurtosis are very sensitive to the 
choice of the zone of the ring. For example, if we were 
to look at the zone kmin = 825, /cmax = 1634 occupied 
by the clusters C~10, we get the skewness and 

kurtosis shown in Fig. 8 (bottom plot). From these two 
plots one can see that the value of the parameter R 70 
is preferred for this analysis. For this scale of filtering the 
kurtosis of the smoothed signal reaches a maximum, while 
for the Gaussian tail of the total signal both the character- 
istics are close to zero both for the whole ring and for the 
particular zone around the CS. In Fig. 9 is shown the ILC 
III signal before (the red dotted line) and after subtraction 
of the smoothed over R c:^ 70 pixels (the black line). One 



sees that the CS is eliminated. 
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Fig. 9.— Top. The ILC III ring b = -57° before demodula- 
tion (red dots) and after (black line). Middle. The ILC III sig- 
nal, smoothed over R = 70 pixels (the corresponding angular scale 
Acf) = 12.3°.) Bottom. The smoothed ILC signal (thick black hue) 
in comparison with signals from Ka (red line), Q (green line) and V 
band (blue line). 



Fig. 10. — The Ka (left) and W band signal from the South hemi- 
sphere. The red circle marks the location of the point source at the 
edge of the CS. 

However, from the middle plot of Fig. 9 we can see that 
the CS, as well as all the large clusters come from the 
smoothed signal, which preserves all the non-Gaussian fea- 
tures mentioned in the previous section including the CS. 
The bottom plot of Fig. 9 show one more significant pecu- 
liarity of the CS zone. Just at the edge of the CS one can 
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see a bright point source. The amphtude of this source is 
about bOcTcmb for the K band and its drops down to Qcicmb 
for V and W bands. In Fig. 10 we show this point source 
for the Ka and W bands. For the ILC III and WCM maps 
the residuals from this point source are associated with 
local maxima at the level of lacmb- 
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5. ARE THERE OTHER "COLD" OR "HOT" SPOTS ? 

This question seems to be one of the most important 
question of the entire analysis of the origin of the CS. Do 
we have the only one single peculiar zone in the CMB 
map, or are there more zones with non- Gaussian prop- 
erties of the signal ? To answer this question, we sug- 
gest using the distribution of the variance of the CMB 
signal versus the latitude, shown in Fig. 5. Let us fo- 
cus attention on the Southern hemisphere and in partic- 
ular on the point of maxima of the variance. We can 
see the following coordinates of the points of maxima: 
= -80°, -70°, -57°, -30°. For these rings we performed 
the same analysis as for the ring at b = —57°. All of 
them show the existence of non-Gaussian features, which 
can be identified by the same method as we have used 
for the b = —57° ring. In Fig. 11 we show the signal 
for the ring —70°. A non-Gaussian hot spot clearly seen 
at A: ~ 850 (the corresponding Galactic coordinates are 
6 = -70°, = 149°). The amplitude of the peak is about 
4.8cro, and is is the member of the cluster with (i+ = 6 for 
Ut = 0. From Fig. 11 one can obtain the optimal size of 
the filter scale R :^ 50. 

Fig. 12 shows the same characteristics of the signal as the 
previous one, but for the b = —30° ring. Once again, the 
characteristic scale of filter is about R = 70 pixels. One 
may continue the search for other rings belonging to the 
Southern hemisphere, just by following the distribution of 
maxima of the variance from Fig. 5. However, the question 
is could we find the same peculiarities of the signal for 
the Northern hemisphere, where a deficit of the variance 
occurs? In Section 3 we already mentioned that the ring 
b = 57° does indeed reveal over-clustering. Let us look 
closely at the skewness and kurtosis for that ring and find 
out the characteristic scale of non-Gaussianity. What is 
interesting is that for the 6 = 57° ring the skewness and 
kurtosis shown in Fig. 13 are very close to the Gaussian 
characteristics, while for the zone kmin = 1520, /cmax = 
1940 with clusters C~15,/cmin = 1520, /cmax = 1731 and 
C+16,/cmin = 1732, /cmax = 1940 wc cau clearly see two 
minima for the demodulated signal, one at = 20 and 
another at = 70. 




Fig. 11. — Top: The skewness and kurtosis for the whole b = —70° 
ring. The definition of the fines is the same as in Fig.9(top). Middle: 
The ILC III ring b = —70° before demodulation (the red dots) and 
after (the black line). Bottom. The ILC III signal, smoothed over 
R = 50 pixels (the corresponding angular scale Acf) = 8.8°.) 
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Fig. 12. The same as Fig.9, but for b = -30° ring. The scale 
of filtration is = 70 pixels. 
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Fig. 13. — Skewness and kurtosis for the b = 57° ring. Top plot 
show these characteristics for the whole ring, bottom plot is for the 
combination of the clustersC~15 + C+16. The color of the lines is 
the same as for Fig. 9. 

However, we stress that for the whole ring at b = 57° 
even without filtering, the skewness and kurtosis are close 
to the GRF, unlike, for example b = —57° ring. Moreover, 
increasing the scale R results in greater departures from 
Gaussianity, as it seen in Fig. 13. This tendency is common 
for b = 72° ring (the point of local minima of the variance), 
as for the b = 78° (the point of local maxima of the vari- 
ance from Fig. 5). Thus, significant non- Gaussianity of 
the ILC III and WCM maps mainly corresponds to the 
South hemisphere, and is associated with large angular 
scales, around 9 — 12°. 

To show that clustering of the extrema of the ILC III 
and WCM signals is a typical feature of the morphology, 
we show in Fig. 14 the the ILC III map seen from the North 
and the South Galactic poles. From the ring with the lat- 
itude b :^ 80° one can easily find a big cluster of maxima 
and one big cluster nearly at the same latitude. For the 
South pole the zone of the CS connect with three big clus- 
ters of maxima. 

6. COMPARISON WITH PECULIAR ZONES DETECTED BY 
WAVELETS 

Detection of peculiarities in the CMB sky is obviously 
one of the major steps in the investigation of the depar- 
ture of the signal from statistical isotropy and homogene- 
ity. As was mentioned in the introduction, this problem 
was discussed by McEwen et al. (2006b), where Spheri- 
cal Mexican Hat Wavelet (SMWH), Elliptical SMHW, and 
Spherical Butterfiy Wavelet (SBW) approaches were ap- 
plied to the detection and location of the position of the 
non-Gaussian spots in the CMB sky. Below we refer to 
these methods collectively as wavelet methods. McEwen 
et al. (2006b) summarize the detection of peculiar zones 
with approximate (estimated) coordinates of these zones. 



We present these data in Table 1 in order to compare these 
results with our results of detection of non-Gaussian zones 
by implementation of cluster analysis. 




Fig. 14. — The image of the sky seen from the North and the South 
galactic poles (top and middle panels). The black circle marks the 
position of the CS. Bottom panel show the South hemisphere of the 
CMB sky after subtraction of the first 20 multipoles of the ILC III 
signal. 
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Table 1 

The peculiar zones of the CMB sky detected by 

WAVELETS (first AND SECOND COLUMNS) FROM McEWEN 

ET AL. (2006b). The second column shows the 

LONGITUDE AND THE LATITUDE OF EACH ZONE. ThE 

third column shows the longitudes of the clusters 
where the wavelet peculiar point is the member. 
The forth column shows the number of extrema in 
the cluster {s) and the length of the cluster {d) . 
The sign means that the wavelet peculiar zone 

IS NOT DETECTED BY CA WITH Z^t = 0. StARS REDIRECT 
THE READER TO FiG. 15. 



Zone 


Location, W 
0, e (deg) 


Location, CA 


S, D 


1 


75, 57 


56, 86 
157, 183 
267, 304 
304, 341 


3, 0.083 

3, 0.072 

4, 0.102 
4, 0.102 


2 


75, 53 


69, 93 
40, 66 

160, 183 
273, 310 


2, 0.067 
4, 0.072 

3, 0.064 

4, 0.103 


3 


323, 56 


304, 339 
278, 304 


4, 0.098 
3, 0.072 


4 


321, 62 


306, 333 
333, 357 
272, 307 


2, 0.075 

3, 0.067 

4, 0.097 


5 


267, 50 


251, 274 

58, 82 
313, 344 


3, 0.064 

3, 0.067 

4, 0.086 


6 


268, 45 


129, 156 


3, 0.075 


7 


213, 40 


30, 58 
233, 258 


4, 0.078 
4, 0.069 


8 


223, 30 


9, 37 
136, 157 


4, 0.078 
4, 0.056 


9 


160, 26 


174, 213 
4, 32 


6, 0.108 
4, 0.078 


10 


94, -28 


35, 62 
142, 235 
333, 352 


3, 0.075 
15, 0.258 

4, 0.056 


11 


81, -34 


141, 181 
187, 230 
261, 294 


7, 0.11 

6, 0.12 

7, 0.092 


12 


118, -42 


8, 30 
333,360+8 


4, 0.061 

5, 0.097 


13 


20, -48 


7, 30 
70, 105 
243, 278 


4, 0.058 
4, 0.097 
4, 0.097 


14 


34, -31 


32, 55 
288, 330 


5, 0.064 

6, 0.117 


15 


230, -68 


220, 253 
19, 52 
63, 117 


2, 0.092 
2, 0.092 
5, 0.15 


16 


204, -56 


198, 216 
148, 188 
223, 270 


2, 0.05* 

4, 0.11 

5, 0.13 


17 


186, -54 


193, 216 
150, 193 
224, 294 


2, 0.063* 
4, 0.12 
6, 0.194 


18 


218, -33 


187, 230 
258, 294 


3, 0.12 
6, 0.10* 






Fig. 15. — The CMB signal for equal latitude rings, marked in 
Table 1 by sign "— " . From the top to the bottom and from the left 
to the right are the zones 6-12,16-18. 

Let US discuss the properties of the signal for Zone 6 of 
Table 1 . The wavelet approach gives us the aziniuthal 
coordinate of the non-Gaussian spot = 268°. From 
the top left plot, we can see that this point corresponds 
to a relatively small negative peak. The same morphol- 
ogy of the peaks is clearly observed at ~ 210°, and 
~ 330°. Both peaks are very similar to the CS zone, but 
with smaller amplitudes. These two peaks probably con- 
tribute to the overall non-Gaussianity of the signal more 
significantly than the 06 = 268° peak (see Table 1). 

The next zone, which was not detected by cluster anal- 
ysis at z/t = 0, is Zone 7 shown on the top right plot. 
According to the wavelet analysis, the azimuthal coordi- 
nate for this zone is 07 = 213°. Once again, here we have 
a negative peak, but the peak at 07 ^ 240° reveals a more 
significant departure from Gaussianity, it being a member 
of the cluster with number of minima N = 4. 

For Zone 8, the wavelets analysis gives us 08 = 223°. 
We show the corresponding ring in Fig. 15 (second from 
the top left plot). The detected zones manifest themselves 
as points of maxima, while the negative peaks detected by 
cluster analysis, listed in Table 1 reveal significant depar- 
ture from Gaussianity. For example, the cluster of minima 
at 140 — 160° can be detected not only for ut = 
threshold, but even for z^t = — 2, as a cluster of maxima 
and minima. Roughly speaking, all the pictures, shown in 
Fig. 15 clearly demonstrate that implementation of cluster 
analysis with different thresholds i^t allow us to detect not 
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only single non-Gaussian zones, but also clusters of these 
zones which have non-local structure. To illustrate this 
tendency, we would like to draw attention to the bottom 
right picture belonging to the Zone 18. According to Ta- 
ble 1, the non-Gaussian zone is located at 0i8 = 218° being 
the member of the cluster and detected both by wavelet 
and by cluster analysis. However, there is another zone 
at (/) c± 80° with a morphology similar to the morphology 
of the CS. Moreover, looking at the shape of the signal at 

> 130 one can see the modulation of the signal by low 
frequency harmonics. This type of non-Gaussianity is an 
argument in favor of hypothesis that low multipoles of the 
CMB signal are highly non-Gaussian. 

7. CONCLUSION 

We have re-examined the properties of the Internal Lin- 
ear combination WMAP CMB map and the co-added 
WCM map by an analysis of the properties of the signal 
in the vicinity of the CS. These two maps of the CMB sig- 
nal display remarkably similar structures on equal latitude 
rings at \b\ > 30°. We have re-examined the properties of 
the CS at the galactic latitude b = —57° and longitude 

1 = 209° and shown that it is associated with the cluster 
with length D ~ 3{D{n)). In addition to the CS, we have 
also found a few more zones of the CMB signal with almost 
the same morphology, at 6 = 57°, 6 = -80°, b = -30°. 

From an analysis of the ILC III map we have shown that 
the shape of the CS is formed primarily by the CMB signal 
localized in multipoles between 10 < ^ < 20 (correspond- 
ing to angular scales about 5 — 10°), in agreement with 
Cruz et al. (2005, 2007) results. Taking into account that 
the same modes lead to a modulation of the whole CMB 
sky, we subtracted these modes from the CMB signal. The 
demodulated CMB signal looks like a random one without 



significant over-clustering. 

We have investigated the asymmetry of the variance for 
iso-latitude rings in respect to the Galactic plane. The 
South hemisphere has excess variance in comparison to 
the North hemisphere. This is why local defects and large 
clusters, including the CS and its associated cluster, are 
mainly concentrated in the Southern hemisphere. 

Taking all these investigations together, we believe that 
the mystery of the WMAP CS directly reflects directly the 
peculiarities of the low-multipole tail of the CMB signal, 
rather than a single local (isolated) defect or a manifesta- 
tion of a globally anisotropic model. This interpretation 
does not preclude the possibility of an exotic origin of the 
CS and related phenomena, but it does specify more pre- 
cisely what properties such explanations must generate. A 
satisfactory model of the CS must explain the entire range 
of its behaviour rather than only one aspect. 

Our flnal remark is related to the deflnition of signifl- 
cance of the CS detection by different methods, based on 
the assumption that Gaussian statistics apply to the ob- 
served CMB sky. Ever since Eriksen et al. (2004), showed 
that the distribution of the power of the CMB across the 
sky is very anisotropic at the scales about 10°, it has 
been clear that Gaussian statistics are no longer a valid 
reference for determining the significance of this feature. 
Our approach to the large-scale angular modulation of the 
CMB is a possible alternative approach to this issue. 
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